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We report a measurement of the rate of prompt diphoton production in pp collisions at -^s = 
1.96 TeV using a data sample of 207 pb~^ collected with the upgraded Collider Detector at Fermilab 
(CDF II). The background from non- prompt sources is determined using a statistical method based 
on differences in the electromagnetic showers. The cross section is measured as a function of the 
diphoton mass, the transverse momentum of the diphoton system, and the azimuthal angle between 
the two photons and is found to be consistent with perturbative QCD predictions. 

PACS numbers: 13.83.Qk, 12.38.Qk 



Diphoton (77) final states are a signature of many in- 
teresting physics processes. For example, at the LHC, 
one of the main discovery channels for the Higgs search 
is the 77 final state 0, Q • An excess of 77 production 
at high invariant mass could be a signature of large ex- 
tra dimensions , and in many theories involving physics 
beyond the standard model, cascade decays of heavy new 
particles generate a 77 signature 0. However, the QCD 



production rate is large compared to most new physics, 
so an understanding of the QCD production mechanism 
is a prerequisite to searching reliably for new physics in 
this channel. In addition, the two-photon final state is 
interesting in its own right. Due to the excellent energy 
resolution of the CDF electromagnetic (EM) calorime- 
ter, the 4- momenta of the two photons in the final state 
can be determined with good precision. This allows, for 
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example, a direct measurement of the transverse momen- 
tum of the 77 system {qt), which is sensitive to initial 
state soft gluon radiation. 

In perturbative Quantum Chromodynamics (pQCD), 
the leading contributions are from quark anti-quark an- 
nihilation {qq 77) and gluon-gluon scattering {gg — > 
77). The latter subprocess involves initial state gluons 
coupling to the final state photons through a quark box; 
thus, this subprocess is suppressed by a factor of with 
respect to the qq subprocess. However, the rate is still 
appreciable in kinematic regions where the gg parton lu- 
minosity is high, such as at low 77 mass. Because the 
probability for a hard parton to fragment to a photon is 
of order aem/cts, processes involving the production of 
one (zero) prompt photons and one (two) photons orig- 
inating from parton fragmentation are also effectively of 
leading order (LO). Next-to- leading order (NLO) contri- 
butions include real and virtual corrections to the above 
subprocesses. 

We have compared our experimental results to three 
predictions : DIPHOX 5], ResBos 0, and PYTHIA 0. 
DIPHOX is a fixed-order QCD calculation that includes 
all of the above subprocesses at NLO (except for gg 77 
which is present only at LO). Recently, NLO corrections 
for gg 77 have been calculated and we have added 
these corrections to the DIPHOX prediction. The Res- 
Bos program includes subprocesses where the two pho- 
tons are produced at the hard-scattering at NLO and 
fragmentation contributions at LO; but also resums the 
effects of initial state soft gluon radiation. This is par- 
ticularly important for examination of the 77 qx distri- 
bution, which is a delta function at LO and divergent as 
(7t — > at NLO, and thus requires a soft gluon resumma- 
tion in order to provide a physical description of the 77 
data in this region. PYTHIA is a parton shower Monte 
Carlo program that contains the above processes at LO. 

At hadron-hadron colliders, it is difficult to measure 
a fully inclusive 77 cross section due to the large back- 
grounds from quarks and gluons fragmenting into neutral 
mesons which carry most of the parent parton's momen- 
tum. Isolation requirements are typically used to reduce 
these backgrounds. In this analysis, the isolation cri- 
terion requires that the transverse energy (Ex) sum in 
a cone of radius R = 0.4(in r] — (j) space) [2j about the 
photon direction, minus the photon energy, be less than 1 
GeV. This isolation requirement reduces the backgrounds 
from neutral mesons decaying into photons and photon 
production from fragmentation sources. The CDF iso- 
lation requirement effectively removes all contributions 
where both photons originate from fragmentation sub- 
processes. However, as will be noted later, some indi- 
cation of single fragmentation subprocesses can still be 
observed in the CDF data. 

The CDF II detector is a magnetic spectrometer which 
is described in detail elsewhere The central detec- 
tor consists of a silicon micro-strip vertex detector in- 



side a cylindrical drift chamber, both of which are im- 
mersed in the 1.4 T magnetic field of a superconducting 
solenoid. Outside the solenoid is the central calorime- 
ter which is divided into an electromagnetic compart- 
ment (CEM) on the inside and hadronic compartment 
(CHA) on the outside. Both calorimeters are segmented 
into towers of granularity Ary x A0 « 0.1 x 0.26. The 
CEM consists of a scintillator-lead calorimeter along with 
an embedded multi-wire proportional chamber (CES) 
located near shower maximum at 6 radiation lengths. 
The CES allows for a position determination of the EM 
shower and for a measurement of the lateral shower 
profile. The average energy resolution of the CEM is 
a{E)/E = 13.5%/ VEsh^ (with E in GeV) and the po- 
sition resolution of the CES is 2 mm for a 50 GeV pho- 
ton. Another important component for this analysis is 
a preshower wire chamber (CPR) mounted between the 
magnet coil and the CEM, at about 1.2/sin0 radiation 
lengths. The CPR detects photon candidates that have 
converted in the magnet coil and other material in the 
inner detector. 

This analysis uses events collected with a trig- 
ger that requires two photon candidates with Et greater 
than 12 GeV each. A requirement of Et greater than 14 
GeV (13 GeV) for the leading (softer) photon candidate 
in the event is imposed in the offline analysis. The mini- 
mum transverse energy requirements for the two photon 
candidates are different in order to avoid the kinematic 
region where the NLO calculation is unstable due to the 
imperfect cancellation of the real and virtual gluon di- 
vergences. 

In identifying photons, we impose fiducial require- 
ments to avoid uninstrumented regions at the edges of 
the CES; as part of this criterion we require the pseu- 
dorapidity of the photon candidate to be in the interval 
1 77 1 < 0.9. The reconstructed z- vertex for the collision 
is required to be less than 60 cm from the center of the 
detector. The ratio of the hadronic energy to EM en- 
ergy (Had/EM) for the photon candidates must be less 
than 0.055-^0.00045 x£:, with E the EM energy in GeV. 
The isolation energy is required to be below 1 GeV. Al- 
though only about 1% of showers from prompt photons 
have more than 1 GeV of additional energy in the iso- 
lation cone, about 15% of the photon showers fail the 
isolation requirement because of additional energy from 
the underlying event |l2j . Photon candidates with any 
tracks {pr above 0.5 GeV) that can be extrapolated to 
them are rejected. The lateral profile of EM showers in 
the CES is compared to the profile of electrons measured 
in a test beam. The definition of the from the compar- 
ison can be found in Ref. . We require the of the 
comparison to be less than 20 in the event selection and 
reject photon candidates with an additional CES cluster 
above 1 GeV ■ The efficiencies for each event selection 
requirement, evaluated using a combination of PYTHIA 
Monte Carlo and data, are listed sequentially in Tabic 
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TABLE I: The selection efficiencies per diphoton event. 

Trigger efficiency 0.951 
Reconstruction efficiency and fiducial 0.423 
Isolation energy in 0.4 cone < 1 GeV 0.727 
No track pointing to the EM cluster 0.699 
No extra CES cluster above 1 GeV 0.899 
CES < 20 0.970 
Had/EM <0.055 + 0.00045x£' 0.976 
I z- vertex I < 60 cm 0.877 



TABLE IL A comparison of the cross section as a function 
of the 77 mass for the data and predictions from DIPHOX, 
ResBos and PYTHIA. 



Combined (etot) 



0.152 



M-y-y 


CDF Data 


DIPHOX 


ResBos 


PYTHIA 


(GeV) 


(pb/GeV) 


(pb/GeV) (pb/GeV) (pb/GeV) 


10-25 


0.03 ± 0.03 ± 0.01 


0.04 


0.01 


0.01 


25-30 


0.44 ± 0.13 ± 0.12 


0.41 


0.31 


0.18 


30-35 


0.61 ± 0.17 ± 0.16 


0.70 


0.65 


0.38 


35-45 


0.46 ± 0.10 ± 0.14 


0.46 


0.43 


0.24 


45-60 


0.16 ± 0.05 ± 0.04 


0.19 


0.16 


0.09 


60-100 


0.01 ± 0.02 ± 0.01 


0.04 


0.04 


0.02 



The trigger efficiency per photon, measured using a sin- 
gle photon trigger, is approximately 80% at 13 GeV and 
rises to greater than 99% for photons above 15 GeV. The 
combined selection efiiciency (etot), including acceptance 
and trigger efficiency, is 15.2% per diphoton event. 

After imposing all of the requirements, 889 two-photon 
candidates remain in our data sample. This sample in- 
cludes background from neutral mesons such as tt" and 
rj that decay to multiple photons. To estimate this back- 
ground, we apply the statistical background subtraction 
method described in which makes use of the differ- 
ences on average between EM showers produced by single 
photons and by the multiple photons produced in neu- 
tral meson decays. The separation between single and 
multiple photon showers relies on the shower shape mea- 
sured by the CES and the preshower conversion pulse 
height measured by the CPR. Since photons from the de- 
cay of neutral mesons with Et above 35 GeV are almost 
coUinear in the lab frame, their shower shape in the CES 
is no longer distinguishable from a single-photon shower. 
To estimate the background contamination in this high 
Et region, the CPR has been utilized. The chance for a 
conversion to take place in the tracking volume or mag- 
net coil (1.1 radiation lengths) and generating a hit in 
the CPR is higher for the multiple photons than for a 
single photon. We use the CES shower shape for photon 
showers with Et < 35 GeV and the CPR pulse height 
for Et > 35 GeV. For each photon candidate, we test 
whether the CES is less than 4 (low Et) or the pho- 
ton canidate produces a pulse height in the CPR greater 
than one minimum ionizing particle (high Et)- There are 
four possibilities for the final state: both candidates pass 
the test, the first candidate passes and the second fails, 
the first fails and the second passes, or both candidates 
fail (the first candidate has the higher Et)- From the 
known efficiencies for photons and background to pass 
the a-nd conversion tests, we can then determine the 
number of true 77 events (as well as the number of 7- 
background, background-7 and background-background 
events). Using the two background techniques discussed, 
we determine that of the 889 candidates, 427 ± 59{stat) 
are real 77 events. From these events, the calculated ac- 



TABLE HI: A comparison of the cross 
the 77 qr for the data and predictions 
and PYTHIA. 



section as a function of 
from DIPHOX, ResBos 



qr 
(GeV) 



CDF Data 
(pb/GeV) 



DIPHOX 

(pb/GeV) 



ResBos 
(pb/GeV) 



PYTHIA 
(pb/GeV) 



0-1 


0.70 ± 0.30 ± 0.14 


-2.45 


0.34 


0.53 


1-2 


1.18 ± 0.43 ± 0.28 


5.59 


0.95 


1.15 


2-4 


0.92 ± 0.35 ± 0.28 


2.06 


1.03 


0.94 


4-8 


0.96 ± 0.23 ± 0.32 


1.17 


0.94 


0.46 


8-12 


0.29 ± 0.21 ± 0.13 


0.44 


0.59 


0.21 


12-16 


0.42 ± 0.14 ± 0.12 


0.24 


0.36 


0.12 


16-24 


0.19 ± 0.09 ± 0.05 


0.13 


0.19 


0.07 


24-32 


0.12 ± 0.06 ± 0.03 


0.09 


0.07 


0.03 


32-40 


0.10 ± 0.05 ± 0.05 


0.06 


0.03 


0.01 



ceptance and the integrated luminosity, we determine the 
diphoton cross sections for several kinematic variables. 
The 77 mass distribution is shown in Fig. ^ along with 
predictions from DIPHOX, ResBos and PYTHIA. The 
qt distribution is shown in Fig. 13 and the A(f> distri- 
bution between the two photons is shown in Fig. The 
vertical error bars on the data indicate the combined sta- 
tistical and systematic uncertainties with the inner tick 
marks indicating the statistical uncertainty alone [lij . 
The PYTHIA predictions have been scaled (factor of 2) 
to the total measured cross section in all the figures. The 
cross sections as a function of these three different vari- 
ables are also tabulated in Tables HTl IIIII and II VI 

It should be noted that the background to the 77 signal 
has been determined independently for each kinematic 
bin as the background fraction can vary with the kine- 
matics. Determining the background on a bin-by-bin ba- 
sis increases the statistical uncertainty but decreases the 
systematic uncertainty. 

The systematic effects include uncertainties on the se- 
lection efficiencies (11%), uncertainties from the back- 
ground subtraction (20-30%) and from the luminosity 
determination (6%) 17] . 

We note some features of the theoretical predictions. 
The gg subprocess provides a significant contribution to 
diphoton production at low mass (^30 GeV/c^); the 
slight wiggle observed in the DIPHOX prediction is due 
to the very rapid falloff (with diphoton mass) of the gg 
subprocess compared to the qq subprocess. The Res- 
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TABLE IV: A comparison of the cross section as a function 
of the 77 for the data and predictions from DIPHOX, 
ResBos and PYTHIA. 



A0-y^ 


CDF Data 


DIPHOX 


ResBos 


PYTHIA 


(tt rad) 


(pb/rad) 


(pb/rad) 


(pb/rad) 


(pb/rad) 


0.0-0.2 


1.06 ± 0.52 ± 0.34 


0.69 


0.01 


0.02 


0.2-0.4 


0.89 ± 0.52 ± 0.32 


0.56 


0.23 


0.09 


0.4-0.6 


0.51 ± 0.63 ± 0.19 


0.71 


0.73 


0.44 


0.6-0.8 


3.34 ± 1.10 ± 1.04 


1.83 


3.08 


1.09 


0.8-1.0 


15.56 ± 2.59 ± 4.70 


23.37 


17.52 


10.68 



Bos qT curve is smooth for the entire range, while the 
DIPHOX curve (as can be seen by the negative value in 
the first bin of Table ITTl|l is unstable at low qT due to the 
singularity noted earlier, and thus is not plotted in Fig. 
Hfor qr <2 GeV/c [13. At the high qr end, DIPHOX 
displays a shoulder, a feature absent in the ResBos pre- 
diction. The ResBos curve lies above the DIPHOX one 
at values of the order of 7r/2 but also lies significantly 
below the DIPHOX curve at small Ai/i. 

The observed differences between the predictions are 
expected. The fragmentation contribution in ResBos is 
effectively at LO. Since fragmentation to a photon is of 
order aem/as, some 2^3 processes such as qg gqj, 
where the quark in the final state fragments to a sec- 
ond photon, are of order a^.^Us and are included in a 
full NLO calculation. These contributions are present in 
DIPHOX, but not in ResBos, which leads to an under- 
estimate of the production rate in the latter at high qx, 
low A(/), and low 77 mass. In particular, the shoulder at 
qx of approximately 30 GeV/c arises from an increase in 
phase space for both the direct and fragmentation sub- 
processes It is instructive to divide the DIPHOX 
predictions into two regions : A<j) > tt/2 and A(/) < 7t/2. 
We do so, and plot the qx prediction for the A<j) < 7r/2 
region in Fig.|21in order to highlight this contribution. It 
is apparent that the bump in the DIPHOX prediction at 
a qx of approximately 30 GeV/c is due to the "turn-on" 
of the A(j> < tt/2 region of phase space. At A(/) values 
above 7r/2, the effects from soft gluon emission (included 
in ResBos but not in DIPHOX) are significant. 

The data are in good agreement with the predictions 
for the mass distribution. At low to moderate qx and Acjj 
greater than 7r/2, where the effect of soft gluon emissions 
are important, the data agree better with ResBos than 
DIPHOX. By contrast, in the regions where the 2— >3 
fragmentation contribution becomes important, large qx, 
A<j) less than 7r/2 and low diphoton mass, the data agree 
better with DIPHOX. 

In this paper, we have presented results for 77 pro- 
duction in pp collisions at a center-of-mass energy of 1.96 
TeV using a data sample twice that previously available. 
Good agreement has been observed with resummed and 
NLO predictions in different regions of phase space. For 
agreement in all areas, however, a resummed full NLO 
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FIG. 1: The 77 mass distribution from the CDF Run II 
data, along with predictions from DIPHOX (solid), ResBos 
(dashed), and PYTHIA (dotted). The PYTHIA predictions 
have been scaled by a factor of 2. The inset shows, on 
a linear scale, the total 77 cross section in DIPHOX with 
(solid) /without (dashed) the gg contribution. 
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FIG. 2: The 77 qx distribution from the CDF Run II 
data, along with predictions from DIPHOX (solid), ResBos 
(dashed), and PYTHIA (dotted). The PYTHIA predictions 
have been scaled by a factor of 2. Also shown, at larger qx, 
are the DIPHOX prediction (dot-dashed) and the CDF Run II 
data (open squares : shifted to the right by 1 GeV for visibil- 
ity) for the configuration where the two photons are required 
to have A(j> < 7r/2. 

calculation will be necessary. 
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FIG. 3: The /S.<j> angle between the two photons from the CDF 
Run II data, along with predictions from DIPHOX (solid), 
ResBos (dashed), and PYTHIA (dotted). The PYTHIA pre- 
dictions have been scaled by a factor of 2. The bump at 
A(/> = 0.3 on the DIPHOX curve is unphysical; it is caused by 
the combination of collinear divergence and the anti-coUinear 
cut [13 ■ 
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